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1 Introduction 

Default swaptions are options on default swaps, hence they are often treated by drawing analogies 
with interest rates swaptions, especially as far as their Black-like market pricing formula is con- 
cerned. Indeed, the most widely discussed model for their valuation, the Log-normal Default Swap 
Market (LDSM) model is similar to the Log-normal Swap Market (LSM) model used in interest 
rates derivatives markets. 

Schonbucher (2000) introduces the notion of survival pricing measures by conditioning on no- 
dcfault up to a given maturity, although in contrast with the interest rate models this measure is not 
equivalent to the risk neutral one. The standard market formula for default swaptions is obtained 
by modeling the default swap spread directly as a geometric brownian motion, as summarized in 
Schonbucher (2004). The convenient pricing measure is termed the survival swap measure whose 
associated "numeraire" is the defaultable annuity which may vanish. 

Jamshidian (2004) partly addresses this problem, presenting a more formalized setup of the 
LDSM model, and generalizes the theory to scmi-martingales driven money market account and 
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conditional survival probabilities. Hull and White (2003) presents practical aspects for the imple- 
mentation of the Black market formula and provides empirical estimates of default swap spread 
volatilities for actively traded A-rated names that ranges from %67 to %130 in the period from 
end 1999 to mid-2002. Brigo (2005) introduces various different candidate formulations by using 
alternative definitions of defaultable forward rates and develops a market model leading to the 
standard Black formula under equivalent pricing measures, showing market implied volatilities of 
the order of %50 . 

The variants of the Black (1976) formulae obtained for the LDSM model, by their inherent 
simplicity are particularly convenient for quoting single default swaptions by selecting an appro- 
priate volatility parameter. However, quoting default swaptions for different sets of maturities 
and strikes or more complex instruments consistently with just the implied volatilities given by 
the formula inversion becomes problematic. One is confronted with the need to develop a fully 
specified dynamical model to impose a structure on the joint dynamics of one-period rates or credit 
default swap spreads as it is done in the interest rates derivatives markets with the LFM, or the 
LSM. 

In particular, the valuation of more exotic instruments like Bermudan default swaptions re- 
quires the use of a model that accurately incorporates the term structure of default swap rates, as 
well as the dynamical deformations and movements of this term structure. Indeed, default swap 
rates are subject to large jumps and possibly stochastic volatility effects. In the interest rates 
derivatives models, these features are incorporated more or less successfully by specifying richer 
joint dynamics of the forward LIBOR rates, leading to an explosion of the number of parameters in 
the models. However, contrary to the interest rate markets with their huge number of caps/floors 
and swaptions, the single-name default swap markets are most famous by the very small number 
of traded instruments, rendering the calibration or estimation of any model with a large number 
of parameters unfeasible. 

An alternative approach, more suitable for the current state of the default swaptions markets 
that is explored in this paper calls for modeling the default intensity instead. This is the approach 
followed in Brigo and Alfonsi (2003) with a stochastic default intensity model where both the 
short rate and the default intensity are driven by shifted square root diffusion processes. Brigo and 
Cousot (2006) examines implied volatilities generated by this two-factor shifted square root model 
and characterizes the qualitative behavior of the implied volatilities with respect to the stochastic 
intensity model parameters. The numerical experiments conducted with stylized parameter values 
suggested that this model might be unable to generate large enough implied volatilities. 

Modeling the intensity process automatically imposes a strong structure on the default swap 
spread joint dynamics across different maturities, simplifying the achievement of consistency across 
instruments. Judicious choices of the intensity process can also incorporate jumps and some 
stochastic volatility effects in default swap spreads, and possibly generate plausible defaultable 
term structure evolutions. 

In this paper, we extend the SSRD model of Brigo and Alfonsi (2003) and Brigo and Cousot 
(2006) by allowing for positive jumps in the process driving the default intensity, consistently 
with empirical evidence. We develop and test a semi-analytical formula for single-name default 
swaptions allowing for fast and accurate pricing. The semi-analytical formula is based on the 
celebrated decomposition due to Jamshidian (1989) for the valuation of options on coupon bonds 
in one-factor affine models. We show with numerical experiments that the model generates plausible 
volatility smiles. Given its relative tractability, the model can thus be calibrated to vanilla default 
swaptions in order to price more exotic products. A first attempt at Bermudan default options 
pricing for example is in Ben Ameur et al (2006), where the basic SSRD model is used. The 
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jump-extended SSRJD introduced here could be used as an improved version. 

2 The SSRJD default intensity model 

Default swap rates time series can hardly be reconciled with a geometric brownian motion. Fur- 
thermore, forward default swaps underlying default swaptions are not traded as such, and hence 
delta-hedging can only be done approximately with spot default swap term structures. For risk 
management and control purposes, it is important to recognize the relation between different de- 
fault swaps and swaptions referencing the same credit name. 

Jumps and stochastic volatility could potentially be introduced in the market model by postu- 
lating more appropriate dynamics for the default swap rate. However, this would quickly destroy 
the main feature of the model: its simplicity. Also, in order to value more exotic options, it becomes 
important to incorporate the whole term structure of the default swap rates as well as postulating 
dynamics that can yield appropriate deformations of this term structure in the future. 

A relatively simple candidate for these tasks is the stochastic default intensity approach we 
adopt in this paper. It models the default intensity process instead of the default swap rate, pro- 
viding a plausible approach to consistently model the default swap rates of different maturities. 
Following this approach, Brigo and Alfonsi (2003) proposes a two-factor shifted square root dif- 
fusion model, where both the short rate and the default intensity are assumed to follow possibly 
correlated shifted square root diffusions. The processes are modeled as a sum of a deterministic 
function and a square root diffusion. 

Comparing numerical examples in Brigo and Cousot (2006) and Brigo (2005) we see that it 
is difficult to produce large enough implied volatilities compared to what is implied from default 
swaptions market data or historical volatilities of default swap spreads. Hereafter we present an 
extension to this model, by allowing for positive jumps in the process driving the default intensity. 

We denote the market filtration by G = (^^ )((>□■) and let Q be a risk- neutral probability 
measure. We follow the intensity based approach to default risk modeling and introduce the 
default time as a totally inaccessible G— stopping time t. We further assume the usual structure 
for G, namely that G = F V H, where F = {J-t){t>o) is the filtration generated by the stochastic 
market variables (interest rates, default intensities, etc) except default events and H = {'Ht){t>o) 
is the filtration generated by the default process: Tit ~ cr u <t). It is also assumed that 

there exists a strictly positive F— adapted process {Xt){t>o) such that the process {Mt)(^t>o) given 
by: 

/•t t'tAr 

(2.1) Mt = l{T<t} - / l{r>s}Asds = l{r<t} - / Asds 

Jo Jo 

is a uniformly integrable G— martingale under Q. The process (At )((>□) is referred to as the G 
marginal intensity of the stopping time r under Q or risk-neutral pre-default intensity. This 
setup is commonly referred to as a doubly stochastic Poisson default process or the Cox process 
framework. In the SSRJD model, the intensity Xt is written as the sum of a positive deterministic 
fimction ip{t) and of a positive stochastic process yt'. 

(2.2) At =yt+^/.(t), t>0, 

where tp is a deterministic function of time, and is integrable on finite intervals. The dynamics of 
{yt){t>o) a-re an example of an Affine Jump Diffusion (AJD) (see Duffie et al. (2000), Duffie et al. 
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(2003) ): 

(2.3) dyt = k{^i - yt)dt + v^tdWt + dJt, 

y(o) = yo, 

with the following condition to ensure the process cannot reach zero: 

(2.4) 2k^i > 1/2. 

(Wt){t>o) is a- Wiener process and {Jt){t>o) is a pure jump process with jumps arrival rate a 
and exponentially distributed jump sizes with mean 7 preserving the attractive feature of positive 
default intensity. In other terms, 

Mt 
i=l 

where (Aft) is a Poisson process with intensity a, the Ys are independent of Mt and also i.i.d. 
exponentially distributed with mean 7. All the parameters yo, k, fj,, v, a,"f are also constrained 
to positive values. Since this model belongs to the tractable AJD class of models, the survival 
probability § has the typical "log-affine" shape before default: 



lit,T) = l{,>t}S(t,T) = l{,>t}Ec 



exp 



Xsds \Tt 



l{r>t}IEc 



exp 



[V^s) + ys]ds \Tt 



(2.5) 

where: 

(2.6) 

(2.7) 



\{r>t}A(t, T) exp ^ i,{^ds - B{t, T)y}j , 



Ait,T) 
B{t,T) 



at,T)C{t,T), 

2(eMT-t) _ 1) 

2/i+ (K + /i)(e^('^-*) - 1)' 



with £,{t, T) and C(^, T) given by: 



(2.8) 
(2.9) 



2/iexp(^(T-i)) 

2h+{K + h){e^^T-t) _ 1) 



2/iexp(^i±ft22(T-t) 



2h+{K + h + 2-f){eMT-t) __ 1) 



and where h = ^/ + 2v'^ . 

Note that when 7 = -^3^, the denominator in the exponent of ((t^T) goes to zero, i.e. ly^ — 
2kj — 27^ = 0, leading to potential numerical instabilities due to division by zero. However, one 
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can check in this case that the base of ({t, T) is then equal to one. Thus for robustness of the 
implementation, it is necessary to set C{t,T) = 1 when 7 = -^^-j^. 

The SSRD model is a diffusion-only restriction of the SSRJD model obtained by setting the 
jump intensity a to zero, also resulting in T) = 1 in the survival probability formula. 

For default swap computations we also make use of the formula for the following transform: 



(2.10) 



exp 



Xsds Xrl^t 



which can be expressed after differentiation as 



(2.11) 
where 



dTS{t,T) =S{t,T) 

dTiit.T) = 
dTB{t,T) = 



2h+{K + h) (eMr-t)_i)^(*' 

4fe^e''(^-*) 

[2h+{K + h) {e^iT-t) - 1)]^' 

-2a7 (e''(^-*) - l) 
2/i + (k + + 27) [e^^iT-t) _ 1) 



Again for the SSRD model, the corresponding formulae can be obtained by simply using the fact 
that C(t, T) = 1 and 9tC(*, T) = 0. 



3 Pricing equations for default swaps and swaptions 
3.1 Credit Default swaps 

In this section, we briefly review default swaps pricing and refer to Brigo and Alfonsi (2003) for 
further details. A (credit) default swap is a financial instrument used by two counterparties to 
buy or sell protection against the default risk of a reference credit name. In a default swap signed 
at time t starting at time Ta with maturity Tf,, the protection buyer pays a periodic fee or spread 
Ra.b(t) at the payment dates Ta+i, . . . , (typically quarterly) as long as the reference entity does 
not default. In case of a default occurring at time t with < r < Tf,, the protection seller 
compensates the protection buyer for his loss given default that we assume to be a known constant 
Lq£). In addition, the protection seller receives from the protection buyer the spread accrued since 
the last payment date before default. In the case where t < Ta, the contract is a forward default 
swap, while if i = we are dealing with a spot default swap. 

Default swaps have been shown in Brigo and Alfonsi (2003) to be relatively insensitive to the 
correlation between brownians driving the intensity and interest rate processes when both are 
modeled as SSRD processes, while Brigo and Cousot (2006) confirms that it is also relatively 
insignificant for default swaptions. Furthermore. Brigo and Cousot (2006) finds that the short 
rate volatility has relatively little impact on the valuation of typically traded default swaptions 
characterized by short maturities, thus concluding that the randomness of the short rate adds little 
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value to stochastic intensity models for default swaptions. Therefore, we assume a deterministic 
term structure of interest rates, and denote the price at time t of the default-free discount factor 
for maturity T or risk-free T-zero coupon bond by D(t, T) = cxp ^— Tgds^ . 

From the perspective of a protection buyer, the value at time t denoted by CDS{t, T, R, Lgd) 
of a default swap with a payment schedule T = {Ta+i, . . . , Tf,}, a spread R and a loss given default 
Lgd is given by the following expression: 



(3.1) 

where 
(3.2) 



CDSit,T,R,LGD)=-l{r>t} 



CaA-t) 



_ n 
RCa.bit) + Lgd / D{t, u)duS{t, u)du 



^ a,D{t,T,)S{t,T,) - I \u^T^,3(^.^^_^)D{t,u)duS{t,u)d^ 



and T/3(t) is the first date in the set {Ta, . . . ,Th} that follows t and ai = Ti — Tj_i is the year 
fraction between Ti_i and Ti. 

Hence, the fair spread Ra,b(t) as long as default has not occurred can be computed as the value 
of R that equates the default swap value to zero: 



(3.3) 



^{T>t}Ra.bit) 



L{r>t}- 



Lgd Jrp'' L){t, u)duS>it,u)du 



3.2 Credit Default swaptions 

A default swaption is an option written on a default swap. In the sequel, we will restrict the 
analysis to European payer default swaptions. A payer default swaption entitles its holder the 
right but not the obligation to become a protection buyer in the underlying default swap at the 
expiration of the option, paying a protection fee equal to the strike spread. Most traded single 
name default swaptions arc canceled (or knocked out) at default of the underlying reference name 
if this occurs before the option's maturity. The maturity of the option will typically be equal to 
the starting date of the underlying default swap Ta- That is, the default swaption holder enters a 
spot default swap if she chooses to exercise the option at maturity. 

For the pricing of a default swaption at a valuation date t, the underlying is thus the Ta 
maturity forward default swap with payment dates Ta+i, . . . ,Tt,. The strike K specified in the 
contract is the periodic fixed rate that is to be paid in exchange for the default protection in 
case of exercise, instead of the fair market spread Ra.b{Ta) that will be available at time Ta only. 
The Tq— defaultable payoff can be valued at time t by taking the risk-neutral expectation of its 
discounted value, where the discounting is done using the default adjusted stochastic discount 
factor D{t, T) = cxp ^— J^{rs + Xs)ds^ as shown in DufRe et al (1996). Hence, the payer default 
swaption can be valued as in Brigo and Alfonsi (2003) and Brigo and Cousot (2006): 



(3.4) 



PSO{t, Ta, T, K) = l[r>t}^Q D{t, Ta)CaATa) {Ra,b{Ta) - K)+ \Tt 



A single dynamics for Ra.b leading to a market formula analogous to the one for interest rate 
swaptions is derived, under different assumptions, in Schonbucher (2004) and Jamshidian (2004). 
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Brigo (2005) derives the same formula under different assumptions and sketches the construction 
of a whole market model for a joint family of default swap rates. Assuming that the default swap 
rate Ra,b follows a geometric brownian motion with volatility aa,b, the above approaches allow to 
price the default swaptions using Black-style formulas. Here, we recall the formula for a payer 
default swaption, with self-evident notation: 

(3.5) PSO{t,T,T,K) = l{,>t}C,.b(t)[i?,,b(t)$(di)-A'$(d2)], 

iogM±(r-t)% 

«1.2 = — 



Oa.h^T ~ t 

When faced with the requirement of marking a default swaption position to market or when hedging 
a book of default swaptions, the need for a different model becomes apparent. Indeed the market 
model requires one to input a volatility parameter. If the model could be trusted as providing 
an appropriate description of the world, this parameter (constant across maturities and strikes) 
could be implied from currently traded options. Recognizing that the model is a rather primitive 
approximation, one would expect to observe different volatility parameters for different strikes and 
maturities resulting in a volatility smile (or skew or smirk) . 

However, for a given underlying reference name, there are often only very few different default 
swaptions traded, and quite often the market is limited to the At-The-Money (ATM) options. 
Deducing patterns in a market model context can then be difficult. On the other hand, the SSRJD 
model can be calibrated to a default swap rates term structure and very few default swaptions, and 
the fitted values of the parameters can be used to value different default swaptions consistently, 
under the condition that the model implies meaningful patterns of implied volatilities. 

To derive a semi-analytical formula for default swaptions in the SSRJD model, we use the 
following equivalent (to (j3.4p ) valuation equation: 

(3.6) PSO{t, Ta, T, K) = D{t, r,)EQ [{CDS{Ta, T, K, Lgd))^ IGt 

Brigo and Alfonsi (2003) proposes a formula for solving this pricing equation in the case of the 
SSRD model. The formula is based on the insightful decomposition of Jamshidian (1989), where 
in a 1-factor yield curve model, an option on a portfolio of cash flows is decomposed in a portfolio 
of options on each cash flow, where the strike for each option is judiciously adjusted. In the next 
section, we prove and extend this formula for the SSRJD model. 



4 Analytical formula for default swaptions pricing 

The derivation of the formula follows three main steps. In proposition 14. 11 we rewrite the pricing 
equation (|3.6p in a suitable form for the application of the decomposition, i.e. as an option on 
an integral of multiples of survival probabilities. Then we use our decomposition in corollary 14. 1[ 
resulting in the appearance of an integral of terms that are akin to options on survival probabilities. 
And lastly, we give an explicit formula for these options in proposition 14.21 Note that proposition 
14.11 and corollarv l4.1l are model-independent in the following sense. Proposition 14.11 holds for any 
nonnegative default intensity process (At) and coroUarv 14. II requires the additional assumption of 
a survival probability function that is monotonic in Aq. We only really use the SSRJD dynamics 
to derive the formula in proposition [521 
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Proposition 4.1. The default swaption price satisfies the following formula: 



(4.1) PSOit,Ta,T,K) = l{,>t}7^(t,Ta)E5 



Lgd - / h{u)Ei{Ta,u)du \ \Tt 



where h is defined as: 

(4.2) h{u)^DiTa,u) [LGD{ru + StM) + K{1 ^ {u ~ Tf,(^^)^i)ru)] , 

with 6t^(u) the Dirac delta function centered at T^. 



Proof. Starting from equation (|3.6[) . we substitute the default swap Tq— value from equations (j3.1 
and (|3.2p resulting in the following formula: 



PSO{t,Ta,r,K) = D{t,Ta)EQ 



7f c,,D{Ta,T,)^Ta,T,) - Lgd D(ra, u)9„S(r„, ) l^t 



We can integrate by parts the last integral of the above expression: 



D{Ta,u)du^{Ta,u)du = D {Ta, u)S{Ta, u) 



S{Ta,u)duD{Ta,u)dv 



Ta 



D{Ta, Thp{Ta,Tb) -I- S{Ta,u)duD{Ta 



u)du. 



For the other integral appearing in the default swaption price, we first decompose it in a sum of 
integrals where the limits of integration are the default swap payment dates: 

''-1 |-T, + i 

D{Ta,u){u-T|3^^.^,)_l)^yS{Ta,u)du^Y^ / D{Ta,u){u ~ Ti)duS{Ta,u)du, 



where we used the fact that for Ti < u < Ti+i, r/3(„)_i = Ti. And we can now integrate by parts 
these integrals: 



D{Ta,u){u-T,)duSiTa,u)du - D{Ta,u){u-T,)S{Ta,u) 



Ti 



Ti 



Ti 



D{Ta, u)§{Ta, u)du 
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Using the fact that Ti — Ti_i = ai, we obtam after summation: 



I To, 



D{Ta,u){u - T^(„)_i)9„§(ra, u)du 



J2 a,D{Ta,T,)S{Ta,T,) 



?— a+1 



D{Ta,u)S{Ta, u)du 



Ta 



§{Ta,u){u- Tg^^)^i)duD{Ta,u)du. 



Note that duD{Ta, u) = —ruD{Ta, u), substitute the expressions obtained for the integrals back in 
the original formula, using 



D{Ta, Tb)S{Ta, Tb) - / DiT,, u)S(Ta, u)STMdu, 

JTa 

and finally, using the formula 

lEQ[l{r>T„}'i^Tj5t] = l{r>t}EQ CXp 



(see for example Bielecki and Rutkowski (2001), Corollary 5.1.1 p. 145), we obtain the result of the 
proposition after rearranging. □ 

Jamshidian (1989) decomposes an option on a portfolio of zero-coupon bonds in a portfolio of 
options on the zero-coupon bonds. The rewriting of the pricing problem as in equation (|4.ip will 
now allow us to achieve a similar result. Indeed, the term J^'' h{u)§{Ta, u)du is akin to a portfolio 
of survival probabilities of infinitely many maturities. We also note that survival probabilities 
satisfy the same formulas as zero-coupon bonds where the default intensity plays the role of the 
short rate. Hence, the expectation in equation (|4.1[) can be seen as a put option on a portfolio of 
zero-coupon bonds (although with infinitely many) where the strike is Lgd and the interest rate 
is given by the default intensity A*. Therefore, it is only natural that we are able to decompose it 
as a portfolio of infinitely many options on survival probabilities. 

Corollary 4.1. Assume the short rate ru is nonnegative and hounded by 100%.' < r„ < 1 \/u G 
[0,Tb] and that the spread payments occur at least once a year^ such that < m — rg(„)_i < 1- If 
the following integral is positive 



(4.3) 



LaoDiTa, u)duS{Ta, u; 0) + /^S(T„ u; 0)7?(r„ u) {l - (u - Tg(„)_i)r„) 



du, 



then the default swaption price is the solution to the following formula: 



(4.4) l[^ytyD{t,Ta) 



h{u)E 



exp 



Xsds {S{Ta,u;y*) - S{Ta,u;yTa))'^\^t 



du, 
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where y* > satisfies: 
(4.5) 



n 



h(u)S{Ta, u; y*)du = L 



GD- 



Otherwise, the default swaption price is simply given by the corresponding forward default swap 
value: 

CDSit,T,K, Lgd)- 



Proof. Recall the definition of h{u): 

h{u) = D{Ta,u) [LGD{ru + ^Ttiu)) + A'(l - (u - T>(„)-i)r„)] . 
Since < r„ < 1 and < u — 7)3(u)-i Ij it follows that 

h{u) > 0, for all u. 

Also note that h{u) is a deterministic function that docs not depend on y, while the survival 
probability S{Ta,u;y) given by equation (|2.5|) is clearly monotonically decreasing in y for all Ta 
and u. Hence, 

rn 

/ h{u)S{Ta,u;y)du 



is a monotonically decreasing function of y. Furthermore, it is easy to see from equation (|2.5[) that 

hm / h{u)S{Ta, u; y)du = < Lgd, 

or just recall that S(Ta,u;y) is a survival probability and y is the initial value of the stochastic 
process driving the default intensity. 

We are interested in finding if there exists y* > satisfying equation (|4.5p . Now, recall that 

h{u) = LGDruD{Ta,u) - K{u - Tp(^^-)_i)ruD{Ta,u) + LGD5Tt,{u)D{Ta,u) + KD{Ta,u) 

and note that (integrating by parts): 

/ \uDiTa,u)S{Ta,u)du^l-D{T,,Tb)S{Ta,Tt)+ f ' D{Ta,u)duHTa,u)du. 
Hence, substituting back in the original integral, we obtain the following: 

LGDD{Ta,u)duS{Ta,u) 

+ KS{Ta, u)D{Ta, u) (l - (u - Tp(u)-i)ru) du. 



h{u)S{Ta,u)du = Lgd + 
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so that: 



lim / h{u)§{Ta,u;y)du 



n 



Lgd+ I LGDD{Ta,u)duSiTa,u;0) 
+ A'§(r„ u; 0)DiTa, u) (l - (u - r^(„)-i)r„) 



du. 



Note that (l — {u — '7)3(„)_i)r„) is nonnegativc since r„ < 1 and {u — T^(„)_i) < 1. And since 
S(Ta,u;0) and D{TaTu) are both nonnegativc being respectively the survival probability and the 
discount factor at time Ta for maturity u, it follows that: 

KS{Ta,u;0)D{Ta,u) (l - (u - T^(„)_iK) > for aU u > T,. 

On the other hand, duS{Ta, u; 0) < since the survival probability is decreasing with maturity. 
We can then consider two cases depending on whether the integral 



^GD 



D{Ta, u)duS{Ta,u; 0) + KS{Ta, u- 0)D{Ta, u) (1 - (u - T0(„)-i)r„) 



du 



is negative or not. 

In the first case, i.e. when the integral is negative: 

hm / h{u)E>{Ta,u;y)du < Lgd, 

and then the equation (|4.5p docs not admit a solution in y. However, in this case the payoff of 
the option is Q — a.s. strictly positive and hence the payoff of the option simplifies to a forward 
default swap payoff. 

In the other case (i.e. when the integral is nonnegativc), by the intermediate value theorem 
the equation (|4.5p admits a unique solution y* by continuity and monotonicity, and we can replace 
Lgd by /j,'' h{u)S{Ta, u; y*)du in (|4.ip . Since Ei{Ta,u;y) is a monotonic function in y, then the 
terms S(Ta, u; y*) — S(Ta, u; yxS) will be all of the same sign for all values of w, and therefore: 



h{u){S{Ta,u;y*) - S{Ta,u;yTj)du 



h{u)(S{Ta,u;y*) - S{Ta,u;yT^))'^du, 



which we can substitute back in the expression (j4.ip for the default swaption value, and use Fubini's 
theorem to change the order of the integrations, resulting in equation (j4.4p . thus completing the 
proof. □ 

Having decomposed the default swaption price in terms of options on survival probabilities, we 
are left with the task of computing these option values. Indeed, to further compute the quantity 
given in equation (14. 4p . recall that: 

A.sds = / ip{s)ds + / Vsds, 
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and that the survival probabUities S(<, T) satisfy equation (|2.5p . Substituting these in formula 
(|4.4p results in the following expression for the default swaption: 



(4.6) 



l{r>t}-D(i, To) cxp J^"" V(s)ds) 



/J" h{u)A{Ta,u)e- '^(^^''^E [exp (- /f" (e-^(^-")y' - e-s(T„,«)y^„)+ |_^^ 



The above expression is analytic up to an integral if we are able to find a formula for the expectation 
involved. We take up that task in the next proposition where: 



^it,T,yt,<;,g) :=E 



exp 



Vsds 



with nonnegative <^ and g. 
Proposition 4.2. 



(4.7) 

where 



'UiT-t,yt,<^,0)-n{T-t,yt,<;,g), 



(4.8) U{T,yo,,,g) = ia^(T)e-^V'(^)^° - - 

Z TT 



°° e^y° [S'cos(W"yo + v<;) + Rsm{Wyo + vc;) 



■dv, 



with 

(4.9) MT) 
(4.10) 



2gh + (2 + g{h ^ k)) {e'^' - 1) 

2/l+(/t + K + £-l^2)(ghT_l) ' 

2 

2/iexp(i^r) 
2h + {h + K + gv^){eJ-''^ - I) 



2h{\ + £17) exp 



(h'-(K+2^Y){l-%(h+K)) 

2(/i-K-27+e(7(/i+K)-i/2)) 



T 



2h{l + g-f) + [h + K + gu^ + 7(2 + g{h - k))] (i 



,hT 



1) 
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and 



R 

S 

U 

w 

E 
F 



{.r + K^)^e'[E cos{H + D arctan(A'/ J)) - F sm{H + D arctan(A7 J))] , 

K'')^e'^[Fi 
5 + ee'^T _^ <^2hT 



{.f + K'^)^e^\Fcos{H + D aYctan{K / J)) + Esm(H + D arctan{K/ J))], 



N 
N ' 



— (x + y cos — -- arctan I — 

= (x^ + y'^)^ sin arctan (^^^^ > 

2fee'''+")¥ [2fe + (fe + K + QU^){e'''^ - 1)] 
N 

2/ie(''+'''¥w2[e^^-l] 



N 



D = 

G = 

H = 

J = 

K = 

N = 

5 = 

£ — 



_ 2'yK — ' 

Q7r[(2 - g(h + K)){h ~ K ~ 27 - g[iy^ - -f{h + k)]) + v'^jh + - -yjh + k)]] 

a7Tt;[(2 - Q{h + k))[v'^ - ^{h + k)] - + - k - 27 - ^([//^ - 7(^1 + «:)])] 



1 + 



(/l - K - 27 - g[l/2 - 7(/i + k)])2 + w2[i,2 _ + ^)]2 

(e'^^ - + K + 27)(1 + ^7) + +7(/i - K))[g{gj + 1) + ^^7]] 



2/1(1 + £17)2 + 2hv^-/^ 

(e''^ - l)u[27fi: + 27^ - i/^] 
2ft(l + £)7)2 + 2/i?;272 ' 

(2/1 + (/i + K + gu^)[e'''^ - l]f + v^Ae'"'^ - 1]', 
2{h — k) — Iz/^fj + g^v^(h + k) + v'^v^{h + k), 
4k — 4k^£( — 2ng'^u'^ — 2v^v'^n, 
—2(h + k) — Av^ g — g^i''^{h — k) — v^u"{h — k). 



Proof. First note the equivalence between the foUowing events: 

{e-'' >e-<^y^}^{yT><;}. 

Hence: 



^{t,T,yt,<,,Q) 



E 



cxp -giyT 
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We define 11 as follows: 



exp ( -gyx - J Vsds j l{yr><;} 



Christensen (2002) derived a formula for 11 that is analytic up to an integral. His formula is also 
reported in Lando (2004) Appendix E. We recall it below: 



n(T, yo, ^, g) = T^i^iT, 2/0, g) 

Z TT 



1 f°°lm[e^'"-ij{T,yo,-g-iv) 



-dv, 



with 



where and satisfy formulae (|4.10p and (|4.9p respectively. The imaginary part appearing 
above admits an explicit expression as given in the statement of the proposition. 
Since the process yt is a homogenous and markovian jump-diffusion 



E 



exp -gyr - I ysds l{.y^><;} |J^t 



exp -gyr-t 



T-t 



ysds l{yr-t>?} 



n(T-t, yt,<,,g)- 



□ 



In summary, if 



LaoDiTa, u)duHTa,u; 0) + X§(r„ u; Q)D{Ta, u) {l - {u - r^3(„)_i)r„) 



du > 0, 



then it is possible to solve for a positive y* satisfying Jrp'' h{u)S{Ta, u; y*)du = Lqd, and such that 
the default swaption price is given by: 



n 



l[r>t}D{t, Ta)e- ^^'^"^ / h{u)A{Ta, u)e- yuy\B{Ta, u))du. 



On the other hand, if 



LcoDiTa, u)duHTa,u; 0) + KSiTa, u- 0)D{Ta, u) (1 - (u - T>(„)-i)r„) 



du < 0, 



the default swaption is so deeply in the money that the probability of it moving out of the money 
is null. Therefore, in this case the default swaption is equivalent to a forward default swap, hence 
it can be valued by computing the price of the equivalent forward default swap. 
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Integral bound: N 


102 


10^ 


10* 


10^ 


10^ 


10^ 


Numerical integral 


-0.75859 


-0.76983 


-0.77173 


-0.77178 


-0.77178 


-0.77178 



Table 1: Numerical approximation of Jp°° ^-^f{v)dv by J" f(v)dv using adaptive Gauss- 

Lobatto quadrature for yo = 0.005, k = 0.196, 0.065, ly = 0.1594, a =^ 0.5, 7 = 0.025, T-t = 1, 
e = S(0,3), <r = 0.0062 



5 Implementation and numerical results 

To implement the formula presented in the previous section, we need to compute the relevant 
integrals numerically. We first focus on deriving a quadrature formula for computing the integral 
appearing in the formula for 11 : 

e'^y [S cos{Wyo + v<;) + R sin{Wyn + v<i)] 

V 

Define f{v) :— S'cos(W^yo + w^) + Rfim{Wyo + ik). In the following lemma^, we confirm that 
our integrand is continuous and bounded on the interval (0, 00) with finite limits on both ends of 
the interval. 

Lemma 5.1. The function ^ J° f{v) is continuous and bounded for v G (0,oo). Moreover, 

pUyo 

(5.2) hm f{v) = 0, 

00 1) 

(5.3) hm f{v) = C, 

ti— »o V 

where C is a constant depending on the model parameters. 

For a visual view of the integrand ^-p-f(v), we plot it for a given set of parameter values in figure 
p]). For the numerical computation of the integral we use the four-point adaptive Gauss-Lobatto 
quadrature with seven point Kronrod refinement provided by Matlab's "quadl" routine based on 
Gander and Gautschi (2000). Numerical convergence can be verified in table ([T]). Experiments 
-not reported here- against a mid-point trapezoidal and Simpson's quadratures confirmed the 
accuracy of the faster and more convenient adaptive Gauss-Lobatto algorithm. For the outer 
integral appearing in the formula (j4.4[) . some experimentation shows that Simpson's rule with at 
worst two quadrature points per quarterly spread payment period is usually enough for convergence 
of the numerical approximation, while using the spread payment dates as the only quadrature points 
in most cases leads to a good accuracy. 

In figure ^ we present some numerical results for payer default swaption prices for different 
strikes, obtained using the quasi-analytic formula developed^. These are for a homogenous non- 
shifted version of the model with constant short rate. The set of parameters used are reported on 
the figure. 



(5.1) 



i 
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Figure 1: Plot of ^f{v) when v e (0,250] for yo = 0.005, k = 0.196, /x = 0.065, = 0.1594, 
a = 0.5, 7 = 0.025, T-t = l 



1 87. 1 5 - 



4-8.27' 
37.74 



Figure 2: Payer default swaption prices (bps) for different strike values (bps) in the SSRJD model 
with no deterministic shifts and parameters: Ta = ly, Tj, = 5y, r = 0.03, Lgd = 0.7, yo = 0.005, 
K ~ 0.229, ^ = 0.0134, V = 0.078, a — 1.5, 7 — 0.0067. The fair value of the underlying forward 
default swap rate is 204 bps. 
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5.1 Consistency with volatility smile 



In this section, we present some numerical results concerning the behavior of the model for some 
parameter values. Our main focus is on the implied volatility smile that can be generated by the 
model. The model potentially allows one to mark-to-market (or rather mark-to-model) non-ATM 
default swaptions that may be present on a trading book. This task cannot be fulfilled with the 
market model unless we use ATM implied volatility to value all options, which should not be 
acceptable from a risk management perspective. On the other hand, our intensity model can be 
calibrated to the default swap term structure and traded ATM default swaptions to price other 
default swaptions more consistently. 

To visualize the implied volatility pattern that can be generated by the model, we present the 
numerical results obtained with three different values of the vector of parameters. The parameters 
values are collected in table 



Reference 


yo 


K 




V 


a 


7 


Modell 


0.0007 


0.4066 


0.0515 


0.1507 


0.5009 


0.0050 


Model2 


1.3E-06 


0.4851 


0.0457 


0.2000 


0.5009 


0.0050 


Models 


0.005 


0.2281 


0.0134 


0.0782 


1.5000 


0.0067 



Table 2 : Three different sets of parameters values of the SSRJD model. 
We plot in figure ^ the CDS term structures generated by the different models. 







-h 


modell 








moclel2 


250 




X 


models 



9 1 O 



Figure 3: CDS term structures for three sets of parameters values of the SSRJD model. The 
default swaps have quarterly spread payments. The short rate is assumed to be constant: r — 0.03. 
The values of the parameters of the intensity process are given in table 0). 

Implied volatility smiles generated from model prices of payer default swaptions with various 
strikes for these models are in figure jl]). Note that the model implies a plausible upward sloping 
volatility smile. It is also noticeable that the presence of a significant jump component can result 
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in higher impUed volatihties than a model with a less important jump component even when the 
last one results in a steeper CDS term structure. 



model 1 
moclel2 
models 



Figure 4: Generated implied volatility smiles for three sets of parameters values of the SSRJD 
model. The short rate is assumed to be constant: r = 0.03. The default swaptions have a maturity 
of one year and can be exercised at maturity into a default swap of a remaining four years to expiry 
with quarterly spread payments. The values of the parameters of the intensity process are given in 
table 0j and the CDS term structures are plotted in figure 



6 Concluding remarks 

The SSRJD model can fit the current default swap term structure while being consistent with some 
dynamic future deformations and implying a volatility smile for default swaptions. The quasi- 
analytic formula presented in this paper permits fast and accurate pricing of default swaptions. 
Hence, the model could be calibrated to the CDS term structure and a few default swaptions. to 
price and hedge other credit derivatives consistently. 
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Notes 



^Usually spreads are paid quarterly. 

^The detailed proof is tedious but straightforward. It is omitted here but is available upon 
request. 

■^Notice that the jumps in the intensity process can only take positive values. If one thinks in 
terms of zero-mean shocks, the long term mean reversion level of the process including jumps is no 
longer the purely diffusive long term mean fi but the larger ^ as summarized in the following 
equivalent way of writing the dynamics of the process yt'. 

dyt = k(/U H yt)dt + Vy^dWt + [dJt - a^dt) 

K 

where the jump increment dJt has been centered by subtracting its mean. This is why, in particular, 
we find a fair value of the underlying forward default swap rate about 204 bps when both the initial 
condition yo and the basic long term mean n are much smaller. 
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